Hybrid squeezing of solitonic resonant radiation in photonic crystal fibers 
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We report on the existence of a novel kind of squeezing in photonic crystal fibers which is concep- 
tually intermediate between the four-wave mixing induced squeezing, in which all the participant 
waves are monochromatic waves, and the self-phase modulation induced squeezing for a single pulse 
in a coherent state. This hybrid squeezing occurs when an arbitrary short sohton emits quasi- 
monochromatic resonant radiation near a zero group velocity dispersion point of the fiber. Photons 
around the resonant frequency become strongly correlated due to the presence of the classical soliton, 
and a reduction of the quantum noise below the shot noise level is predicted. 



Introduction - In quantum mechanics the variance of 
two conjugate observables must respect a minimum value 
set by the Heisenberg uncertainty principle, known as 
standard quantum limit or shot-noise level. Squeezing 
is the result of strong correlations, generally induced 
by nonlinear optical processes, such that the variance 
of a given observable is smaller than the shot noise 
Squeezed light constitutes an important tool for high- 
precision measurements (2j], as well as for applications 
in quantum cryptography and quantum information pro- 
cessing (see Ref. [3|] and references therein). Two quali- 
tatively different kinds of squeezing are currently known 
to take place in optical fibers. The first is due to the pro- 
cess of four- wave mixing (FWM) Q , in which an intense 
continuous wave (CW) pump field generates signal and 
idler waves at symmetric frequencies around the pump. 
Energy conservation leads to a strict photon number cor- 
relation between the signal and idler twin modes, result- 
ing in squeezing. The second kind of squeezing is due to 
self-phase modulation (SPM) [H'l, in which the spectral 
components of an intense short pulse (e.g. a soliton Q) 
in a coherent state become correlated. This is different 
from the FWM-induced squeezing, in that the input field 
is pulsed and is itself squeezed as opposed to squeezing a 
vacuum state 0, S] ■ 

In this paper we theoretically demonstrate that a third 
kind of squeezing is possible in optical fibers. We in- 
vestigate the quasi-monochromatic dispersive radiation 
emitted by optical solitons subjected to perturbations. 
Evidence of photon correlation is observed, manifesting 
itself in a noise spectrum of the emitted radiation lower 
than the standard quantum limit. 

Resonant radiation amplitude with Green's function 
method - Our starting point is the generalized nonlinear 
Schrodinger equation (GNLSE) in the absence of Raman 
effect: id;,A + D{idt)A+\A\'^A = 0. Here z and t are the 
dimensionless time and space coordinates, scaled with the 
input pulse duration Iq and the second-order dispersion 
length Ld2 = *o/l/^2(wo)|, respectively. A{z,t) is the di- 
mensionless electric field envelope in units of y/Po, with 
Po = 1/32(^0)1/(7*0)1 1^2 is the group velocity dispersion 
(GVD) coefficient at the central frequency uto of the input 



pulse, and 7 is the nonlinear coefficient of the fiber [9]. 
The linear dispersion operator in the GNLSE is given by 
D{idt) = is9?+E™>3 am[«9t]™, where s = +1 (s = -1) 
for anomalous (normal) GVD, am = /3m/[m! |/32|io"~^]j 
where /3m is the m-th order GVD coefficient. When 
no higher-order dispersion (HOD) coefficients are taken 
into account {am = 0,Vm > 3), the GNLSE takes the 
standard form of the conventional NLSE with the follow- 
ing fundamental soliton solution A(z,t) — F{t) eyij>{iqz), 
with F{t) = y/2qsech{y/2qt) , where q is the soliton 
wavenumber Q. In absence of HOD, the fundamental 
soliton is stable and invariant along propagation. When 
the NLSE is perturbed by HOD terms, solitons propa- 
gating according to the GNLSE emit a special kind of 
quasi-monochromatic dispersive radiation, called reso- 
nant (or Cherenkov) radiation (RR), provided that phase- 
matching between soliton wavenumber and the fiber dis- 
persion is met [13, [lH ■ 

We now proceed to calculate the amplitude of RR 
emitted by a soliton in the presence of HOD by using 
a Green's function method, a procedure that is intro- 
duced for the first time in this paper. This will pro- 
vide an analytical expression for the spatiotemporal evo- 
lution of creation and annihilation operators of the RR 
around the resonant frequency. These will be used to 
calculate explicitly the quadrature variance, in order to 
investigate quantitatively the squeezing effect along the 
fiber. We make the following ansatz in the GNLSE: 
A ~ [F{t) -I- 5(2, *)] exp(jgz), where g {g'^) is the quan- 
tum annihilation (creation) operator associated to the 
radiation amplitude, which is much smaller than the 
soliton amplitude. By keeping only the linear terms in 
g = [g.g^)^ we obtain: 



P{^^t) 



g^U{F)~Mg, (1) 



where a — diag(l,— 1) is the third Pauli matrix. The 
soliton potential matrix M{t) and the HOD perturbation 
operator P(idt) are defined by (from now on an overbar 
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indicates complex conjugation): 



M{t) = 



2\F{t)\' 



P{idt) = 







P{idt 


(2) 

where P{idt) = [D2{idt) - D{idt)], = \sdl U{F) = 
iP{idt)F,P{-idt)Ff = P{idt)F, with F = iF,Ff. 
Hermitian matrix M is the potential felt by radiation g 
due to the soliton, while vector U{F) is the radiation's 
source term due to HOD perturbations. Note that be- 
cause M is non-diagonal, field g will be coupled to field 
g^ for z 7^ 0, even if there is no initial coupling at z = 0. 
This feature will eventually be responsible for the squeez- 
ing of RR, which is thus indirectly induced by the soliton 
itself. 

We now solve Eq. ([T|) by first solving the following 
equation for the Green's function Q: 

[id, - q + Dildt)] g{z, t-t') = Ait - t'), (3) 

where A{t) is the Dirac delta function. Once Q is known, 
the convolution theorem allows one to find the general 
solution of Eq. (P) : 

g{z,t) ^ go{t)+ Q{z,t-t') \u{F{t'))~ M{t')g{z,t')\ dt' , 

J —oo 

(4) 

where go is the quantum radiation field at z = (corre- 
sponding to vacuum fluctuations), and we have defined 
the 2x2 matrix g = diag(t/,^). Note that Eq. g]) is 
an integral equation: the unknown radiation field g{z, t) 
appears on both sides. To solve Eq. (j4]), one needs first 
to find Green's function Q by solving Eq. ([3]), and then 
by using an iterative perturbative procedure it will be 
possible to find g{z,t) with an arbitrary accuracy. 

It is possible to rigorously show by using the residue 
theorem of complex analysis that the z-dependent 
Green's function Q is given by 



g{z,t-t') = 



2tt 



+00 g— itjft- 



Jl-q+D{io)]z 



-q + D{uj) 



duj, 
(5) 

where P indicates the Cauchy principal value of the in- 
tegral, and the integrand is an analytic function. An 
analysis of the pole structure indicates that if the soliton 
emits RR due to perturbations, the resonant frequencies 
will be solutions of the equation —q + D{uj) = 0, where 
Dioj) is the Fourier transform of the dispersion opera- 
tor D(idt). We assume in the following (without loss of 
generality) that there will be only one real root of this 
polynomial equation, namely uj = Sq. This occurs easily 
in the context of the third-order dispersion induced RR 
emitted by PCFs llj. Moreover, one should note that 
the integrand in Eq. ([5]) contains an exponential with ar- 
gument i[—q-\-D(uj)\z^ which can be, generally speaking, 
a complicated function of the complex frequency lj. It is 
therefore not possible to solve the integral of Eq. ^ ex- 
plicitly, unless a specific integration path satisfying Jor- 
dan's lemma [l6| is found. However, such a path can be 



found once an additional (but not excessively restrictive) 
assumption is made, i.e. that the function [—q + D{(jj)] 
can be expanded in a Taylor series around the only res- 
onance 



-q^D{uj)] 



[oj ~ So] =Vg^[uj~ So], 



(6) 

where Vg{So) is related to the group velocity at the RR 
frequency. The explicit form oi Q — diag(^,^) deduced 
from Eqs. ([5][6l) is given by 



g{z, t-t')= ivgae'^'°^''-'^,{vg,t' - t), 



(7) 



where f}z(vg,t' - t) = dia%{riz{vg,t' - t),riz{vg,t' - t)), 



and riz{vg,t' — t) = [sgn(i' — t) — sgn(t' — t 



-iz)]/2 



is a useful 'truncation' function that appears often in 
the calculations, and physically provides the position 
of the RR emitted by the soliton in the time do- 
main. We are now in the position to use the convolu- 
tion theorem Eq. (j4]) and evaluate the radiation field 
g{z,t) perturbatively. Let us define a new matrix ker- 
nel K{z;t,t') = —Q(z,t)M{t'), and a 'constant' vector 
Go{z,t) = giz,t- t')U{F{t'))dt\ which physically 
represents the classical background field associated with 
the soliton tail, source of the RR classical parametric 
growth. With these new definitions we can write Eq. ^ 
in the following symbolic way: 



ga + Ga + Ki 



(8) 



where the symbol ® denotes the operation of convolution 
with respect to the time variable. Using this symbolic 
notation, one can rearrange the terms of Eq. ^ to obtain 
a formal solution for the amplitude field in terms of a 
Lippmann-Schwinger equation 15 1 
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-(50 + Go). 



(9) 



We can expand the operator [1 — K(^$] ^ in a series, 
[1 - = l + k(g)+K(g)K(g)+.... The 0-th or- 

der in the perturbation theory is given by tf^^Zjt) — 
go{t) + Go{z,t). At this level the interaction between 
the solitonic potential M and the radiation is totally ne- 
glected, the coupling between fields g and in Eqs. (fT2l 
1131) is not present and squeezing is not possible. Note 
that, however, the 0-th order gives reasonable approxi- 
mate results if one is only interested in the calculation of 
the classical RR amplitude 13, HI- To the next order in 
the Born series, one obtains 

= [l + K(E)]{go + GQ) =gQ + k®go + [Go + K^Ga], 

(10) 

which already takes into account the potential matrix M 
hidden in the kernel operator K. The term in square 
brackets in Eq. ((TU)) represents the influence of HOD 
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acting on the sohton, and is a purely 'classical' quan- 
tity. In this paper the main focus is on the quantum 
noise properties of the RR. Thus we decouple the classi- 
cal background soliton tail Gq from the equations, since 
this does not affect the evolution of the quantum fluctu- 
ations around the resonant frequency. 

Now we can specify the initial quantum vacuum field at 
the resonant frequency So (the frequency detuning from 
the soliton carrier frequency): go{t) = (ao(t), aj(i))"'", 
where aj and do are, respectively, the i-dependent cre- 
ation and destruction operators at the input of the fiber 
(z = 0), which satisfy the shot noise commutation rela- 
tion [ao(i), aj(t')] = A{t — t'). In quantum optics this 
vacuum fluctuation field is everywhere and always ex- 
ists even in the absence of classical signals. Due to Eq. 
()10|) . vacuum noise fluctuations interact with the spec- 
trally separated soliton, and acquire a z-evolution given 



J 



hygi^){z,t)^{a{z,t),a\z,t))'^. 

From the above definition of the kernel operator K we 
have: 



' go 



+00 



g{z,t-t')M{t')goit')dt'. (11) 



By using Eq. ([2]) and Eq. ([7]) the above con- 
volution can be expressed as follows: K ® go = 



dt' . It can be shown 



~Vg{H,H)\ where H 
t) [2\Fit')\^ao{t') + F^it')al{t') 
that for an arbitrary function ip one has Vz i^g j t' ~ 
t)ip{f)dt' ^ J^''"^'" ip{f)dt'. 

Squeezing of resonant radiation ~ By substituting 
Eq.([ni) into Eq. ^ we have: 



d{z,t) = do{t) — ivge 



-iSot 



^{z,t) = alit) + tvge'^°' 



t — z/v 



t-Z/Vg 



[2\Fit')\^ao{t') + F^{t')alit')} dt' 
[F'it')aoit') + 2\F{t')\'alit')} dt' 



(12) 
(13) 



Equations (|12m3l) arc the main results of this paper. 
They show that creation and annihilation operators as- 
sociated with the RR photons become correlated via the 
classical soliton F{t), even if photons are uncorrelated at 
z = 0. This correlation leads to a novel kind of squeezing 
(which we call hybrid squeezing), that is conceptually an 
intermediate case between the SPM-induced squeezingof 
pulses Q and the FWM- induced squeezing of CWs 3j 
and is reported for the first time in this paper. 



In order to investigate the properties of the hybrid 
squeezing we construct the two-time radiation noise func- 
tion V{z,t,t') of the quadrature operator X{z,t) — 
a(z,t)e~** -|- {z,t)e^'^ , where $ is a local oscilla- 
tor phase [j, ||: V{z,t,t') = {0\X{z,t)X{z,t')\0) - 
{0\X{z,t)\0){0\X{z,t')\0), where |0) denotes the vacuum 
state. The full temporal dependence of V must be re- 
tained, since the temporal profile of the soliton breaks the 
translational invariance of V along t, analogously to what 
happens in SPM-induced squeezing J8j]j but unlike the 
case of FWM-induced squeezing Function V{z,t,t') 
can be written in a semi-analytical (but somewhat cum- 
bersome) form, which we do not show here. 

The radiation noise spectrum S{z, lu) is a real function. 



and can be calculated by using the standard formula 

S{z,iu) EE - / dt dt' V{z,t,t')e"^^*-' \ 

(14) 

where T is the integration time of the spectrum analyzer 
. Condition 5* = 1 corresponds to shot noise level, and 
S{z = 0) = 1 at the beginning of propagation. For suffi- 
ciently long integration times (T ^ to), expression (fT4)) 
turns out to be independent of T, since V{z,t,t') is a 
localized function in both t and t' due to the soliton lo- 
calization itself. Squeezing of quantum noise is achieved 
in those spectral regions for which S{uj) < 1. We now 
restrict our attention to the resonant frequency a; = ^o- 
In Fig. dja), S{uj = Jq) as a function of z is plotted 
with fixed values of $ and as, for two different values 
of the soliton amplitude q. This demonstrates the exis- 
tence of an optimal value of the fiber length for which 
maximum squeezing of RR is achieved. For large val- 
ues of z, squeezing is lost, so that one would need rela- 
tively short pieces of fiber to observe the described phe- 
nomenon. This is also advantageous in practice, since in 
this case the Raman effect, which manifests itself through 
the Raman self- frequency shift of solitons ^4i| , does not 
have time to take place and can be neglected, as we also 
do here in the equations. Of course in a more precise 
treatment the effect of the Raman nonlinearity on the 
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FIG. 1: (Color online) Squeezing spectrum in various condi- 
tions and parameters, (a) Noise spectrum ^(^o) versus prop- 
agation distance z for fixed value of the local oscillator phase 
"1? = 1.8. (b) Noise spectrum Si&o) versus $ for fixed value 
of propagation distance z = 0.28. The green dotted and red 
solid curves in (a,b) correspond, respectively, to the soliton 
wavenumber q = 0.5 and q = 0.25. The black strait lines in 
(a,b) represent the shot noise level, (c) Contour plot of {S - 
1) showing the z;-evolution of the noise spectrum at various 
noise frequency lj. Parameters used in (c): $ — 1.8; q — 
0.25. Other parameters used in (a,b,c): as — -0.2; all other 
HOD coefficients vanish; integration time T = 10; with these 
parameters the resonant radiation frequency So ~ -2.7. (d) 
GVD of the PCF used in this work. At the soliton wave- 
length (As = 1.2 fim), as ~ -0.2. 



correlations must be considered. This analysis will be 
reported elsewhere. 

Fig. [TJb) shows the evolution of S{6q) versus the local 
phase oscillator $, for fixed values of z and a^, plotted for 
two different values of q. In this case, due to the structure 
of function V{z, t, t'), there is a periodic dependency of S 
on $, and maximum squeezing occurs periodically in $ 
with a period equal to tt. In most cases, in our numerical 
simulations we managed to obtain states with squeezing 
parameters S{So) down to around 0.9, which should be 
easily detectable in possible future experiments. 

It is also interesting to explore what happens to the 
radiation spectrum around uj = Sq. In Fig. [Ijc) the 
contour plot of S{uj) — 1 versus the propagation distance 
z is shown. For as = —0.2, the solitonic source of the 
radiation is located at a; — Jo = 2.7, exactly at the right 
edge of the spectral window. In the plot, the blue colour 
indicates those regions for which squeezing is possible 
[S{u}) < 1], while red indicates anti-squeezing [S{lu) > 
1]. Such figure also shows that, for each given value of 
Lu, there is an optimal value of z that shows maximal 
squeezing, consistent with Fig. [TJa). Moreover, it is 
evident from Fig. [Tfc) that one has a more pronounced 
squeezing effect in close proximity of the spectral body 
of the soliton. 

Proposed experimental verification - Here we propose 




FIG. 2: (Color online) Schematic of the proposed experimen- 
tal verification. A pump laser provides optical pulses that 
propagate through the PCF as solitons. The generated reso- 
nant radiation (RR) is isolated with a spectral filter (SF) and 
coupled into an optical cavity. A photodiode at one cavity 
output in combination with a spectrum analyzer can be used 
to measure the noise spectrum around the RR frequency Sq. 
Tuning the cavity provides access to the quantum fluctuations 
at any field quadrature. 



an experimental way to investigate the squeezing of RR 
in PCF. As a source, solid-core PCF is advantageous, 
since its GVD can be engineered almost at will by choos- 
ing an appropriate hole spacing (pitch) A and hole size 
d in the fabrication process. This offers great flexibility 
in where the RR can be generated and the fiber struc- 
ture proposed here is merely one of many possible exam- 
ples. In Fig. md) we show the GVD and the refractive- 
index profile of a solid-core PCF with A = 1.1 fim and 
d — 0.8 /im. This fiber exhibits a third-order dispersion 
coefficient as = —0.2 at the input soliton wavelength 
Xs = 1.2 /zm, values that we have used throughout this 
paper. The resulting RR is generated in the normal dis- 
persion regime at Ao — 1.38 /im. To measure the gener- 
ated squeezed vacuum at frequencies around Aq one can 
employ usual homodyne techniques, in which the RR is 
mixed with a strong local oscillator (LO) field that shares 
the same mode characteristics as the RR. Alternatively, 
this mode-matching condition can also be met by em- 
ploying a self-homodyne method based on the disper- 
sive action of optical cavities near resonance [13, EH . In 
this approach, the carrier of the RR essentially acts as 
the LO and is interfered with its sideband frequencies 
via an empty optical cavity. The signal reflected by the 
cavity is sent to a photodiode whose photocurrent can 
be studied using a standard spectrum analyzer (Fig. 
Note that the analysis frequency corresponds to the "fre- 
quency detuning" shown in Fig. [IJc). Quantum fluctua- 
tions at any quadrature are simply accessed by scanning 
the cavity's detuning with respect to Sq. This causes the 
noise ellipse to rotate relative to the mean field and thus 
corresponds to changing the LO phase in a traditional 
homodyne setup. The use of a properly designed PCF 
could also help to suppress the guided acoustic wave Bril- 
louin scattering (GAWBS, see Ref. [l^l), which affects 
the squeezing performance and is neglected here for sim- 
plicity. A more precise analysis will include the effect of 
the Raman and Brillouin scattering in the formalism. 

Conclusions - By using an original method based on 
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Green's function, we have demonstrated the existence of 
a novel kind of squeezing, which occurs in correspondence 
of the resonant quasi-monochromatic radiation emitted 
by a soliton near the zero-GVD point of a PCF. Correla- 
tion between photons occurs at the resonant frequencies 
and is induced by the presence of the effective potential 
generated by the soliton. An experimental configuration 
for observing the presented phenomenon has been pro- 
posed. Important applications of this new kind of squeez- 
ing include quantum-enhanced measurements, quantum- 
information processing and quantum cryptography. 
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